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Abstract 

We compute the order nonperturbative contributions to the inclusive 

differential B — > X c lv decay rate. They are parametrized by the expectation 
values of two local and four nonlocal dimension-six operators. We use our 
results to estimate part of the theoretical uncertainties in the extraction of 
matrix elements A and Ai from the lepton spectrum in the inclusive semilep- 
tonic B decay and find them to be very large. We also compute the 
corrections to the moments of the hadronic invariant mass spectrum in this 
decay, and combine them with the extracted values of A and Ai to put an 
upper bound on the branching fraction Br(B — > D**iv). 



1 



I. INTRODUCTION 



Over the last few years there has been much progress in our understanding of the inclu- 
sive decays of hadrons containing a single heavy quark. Combining heavy quark effective 
theory (HQET), with the operator product expansion (OPE), enabled one to show that the 
spectator model decay rate for B — > XJLv is the leading term in a well-defined expansion 
controlled by the small parameter Kqcd/^q, where ttiq is the heavy quark mass fflj. Non- 
perturbative corrections to this leading approximation are suppressed by two powers of m,Q, 
and are parametrized by the matrix elements 

Ai = (Hooiv) | h v (iD±) 2 h v \ H^v)) , (1) 

and 

A 2 = (ffoo(v) I K 9 - a, u h v | Hjy)) , (2) 

where h v is the quark field in the heavy quark effective theory. | H^v)) is the pseudoscalar 
(dp = 3) or vector (dy = —1) heavy meson state in the infinite quark mass limit |§-fl) 
with normalization (Hoo(v) | H^v')) = (2-k) 3 v°5^ 3 \p — p'). The scale dependent matrix 
element A 2 can be obtained from the measured B* — B mass splitting, A2(m&) ~ 0.12 GeV 2 . 

The determination of quantities like Ai and the b and c quark pole masses from exper- 
iment is complicated by the presence of ultraviolet renormalons 0. If the renormalons are 
present, the values of an HQET matrix element extracted from two different observables 
at a given order in a s may differ by an amount of the order of the matrix element itself 
0, which prevents one from using the measured value of one observable to improve the 
prediction for another. Whether this is the case can be established by expressing the un- 
known HQET matrix element in terms of the first observable and substituting this into the 



*In the "large /3q" approximation Ai does not have a renormalon ambiguity in continuum regu- 
larizations || but this is likely to be an artifact of this approximation. 
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theoretical formula for the second. Only if the resulting expression has a reasonably well 
convergent expansion in powers of a s , it makes sense to use the value of the HQET matrix 
element extracted from the first observable to predict the value of the second. In practice, 
one knows only a few terms in the perturbative expansion, and it is hard to assess how well 
the series converges. 

Recently Ai and the difference between the meson masses and the pole quark masses, 
A, have been extracted from the measured inclusive lepton spectrum in semileptonic B 
decays §: Ai = -0.19 ± 0.10GeV 2 ,A = 0.39 ± O.llGeV. The quoted uncertainties are 
the statistical errors only. There are reasons to think that systematic experimental errors 
are not very large. The major theoretical uncertainties come from order perturbative 
corrections, the assumption of quark-hadron duality, and the higher orders in the heavy 
quark expansion. For a very similar analysis see [Tj|. An independent constraint on A and 



Ai can be obtained from the inclusive hadron spectrum in B decays ||. 

Here we compute the terms of order 1 / mf in the heavy quark expansion of the differential 
decay rate B — > X c £u and use the results of our calculation to estimate part of the theoretical 
uncertainties in the determination of A and Ai from inclusive B decays. There are two sources 
of \jm\ corrections. First, the OPE has to be extended to include the local dimension-six 
operators. Second, the lower order corrections calculated in Refs. [@-§j] are expressed in 
terms of the expectation values of dimension-five operators between the physical B states, 
rather than between the states of the effective theory in the limit — > oo. Therefore 
they depend on beyond leading order. In Sect. II we compute the contributions from 
the local dimension-six operators to both the charged lepton spectrum and the hadronic 
spectrum, which are experimentally accessible quantities. The mass dependence of the states 
is discussed in Sect. III. The complete l/mf corrections are parametrized by the expectation 
values of two local and four nonlocal operators. In Sect. IV we investigate the influence of 
1 / mf corrections on the extraction of A and Ai from both leptonic and hadronic spectra in 
B decays. We also obtain an upper bound on the branching fraction Br(B — > D**£v). Our 
conclusions are presented in Sect. V. The Appendix contains the derivation of the meson 



mass formulas to order AQ CD /rriQ. 

II. LOCAL DIMENSION-SIX OPERATORS 

The effective Hamiltonian density responsible for b — > civ decays is 

H w = —Vcb—j^-J^J^j (3) 

where J M = cl7 m &l is the left-handed quark current, and Jjf = IlI^^l is the left-handed 
lepton current. The differential decay rate is determined by the hadronic tensor 

W^ = (2rr) 3 J2S i ( P B-q-px c )(B(v)\J^\X c )(X c \J^\B(v)), (4) 

x c 

which can be expanded in terms of five form factors: 

= -g^Wi + v»v u W 2 - ie^ aP v a q p W 3 + q^q u W 4 + ( 9 V + q u v")W 5 . (5) 

Then the differential semileptonic decay rate is given by 

dT = 96To (6) 
dq 2 dE t dE v m\ [ ' 

x (Wx q 2 + W 2 (2E t E v - ^ + W 3 q 2 (E t - E v )^j 9(E e ) 9{E V ) 6(q 2 ) 6(4E e E u - q 2 ). 
Here To is the spectator model total decay rate in the limit of zero charm mass 

r. = W^Gl^, (7) 

and we have neglected the lepton mass. 
We define the current correlator T^ v by 

= -i J d 4 xe- iq - x (B(v)\T [r\x)J»(Q)] \B{v)} (8) 
= -g^T x + i;VT 2 - ie^VaqpTs + q^q v T 4 + {q»v u + q u v^)T 5 . 

One can easily see that W,- t = — ^ImTj. Away from the physical cut T^ u can be computed 
using the OPE [|TJ. Then analyticity arguments show that the smeared differential decay 






c c 
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FIG. 1. (a) The relevant term in the operator product expansion. Wavy lines denote the 
insertions of left-handed currents, (b) does not contribute to b — > c decay. 

rate is correctly reproduced by the OPE calculation, provided the width of the smearing 
function is large enough. 

The only diagram which has a discontinuity across the physical cut is shown in Fig. la. 
The corresponding contribution to the time-ordered product is 

h v Pt , -fPiJb (9) 

= —67 P L (m b f - $ + t]p+m c ))^ * tPlO, 

A o n=0 \ ^0 / 

where Pl = |(1 — 75) is the left-handed projector, A = {m h v — q) 2 — m 2 c + iO, D is the 
covariant derivative, and we used D^D V — D U D^ = igG^. The field b(x) in eq. (H) is related 
to the normal QCD field by 6qcd(^) — e~ mbV ' x b(x). There are other contributions in the 
OPE of two currents, e.g., the one in Fig. lb. However these operators do not contribute to 
the decay rate once sandwiched between the S-meson states. For the diagram in Fig. lb 
this is ensured by m c being much larger than the available energy in the "brown muck," 
which is of order IVqcd- 

Our calculation of the form factors Tj follows the method of Ref. 0. We expand eq. (0) 
to third order in D. The term with no derivatives is proportional to the conserved current 
tr/^b, and thus its diagonal matrix elements can be evaluated exactly in full QCD. All other 
contributions we express in terms of the field h v in the effective theory and reexpand the 



resulting expressions in powers of l/m&. Therefore we need the expression for b(x) in terms 
of h v (x) only to order 1/?7?|: 



ip ± , (v-D)p ± pi 

where D± = D — v(v ■ D). We choose to work with Foldy-Wouthuysen-type fields, because 
this ensures that they satisfy the usual equal-time commutation relations [|l(| . 



= 1 + 2^ + ^p-S + - K(X) ' <10) 



To evaluate the expectation values of the heavy quark bilinears we need the equations 
of motion in the effective theory to order 1/t?t| |jT0jJTT[| : 

iv ■ Dh v = (±-Pl - ■ D)P ± + ^ (Pl(v ■ D) + (v ■ D)Pl) + ■ ■ •) h v . (11) 

By virtue of eq. (|ll|) there are no nonperturbative corrections to the form factors Tj at order 
l/ra& 0. The contributions at order are expressed in terms of the matrix elements 

(B(v)\h v (iD ± fh v \B(v)), (12) 
l(B(v)\h v l <v K\B{v)). 

Our calculation of these contributions reproduces the results in Ref. ||. The states in the 
matrix elements eqs. ( |T2"| ) have an implicit dependence on m^. At order 1/m^ this dependence 
can be neglected, in which case these matrix elements may be replaced by Ai 2 defined in 
eqs. (0) and (fj). If the form factors are to be calculated to order this replacement 

is no longer valid. An expression for the matrix elements eq. (|T2"|) in terms of A12 and the 
expectation values of nonlocal operators is given in Sect. III. 

The l/m| contributions to the form factors from local operators can be parametrized 
by two matrix elements, p\ and p 2 |12 Pi They are defined as 



(Jl^(v)\h v {iD a ){iD^(iDp)h v \H^{v)) = \pi {g aP - v a v p ) v^, (13) 
(H ca (v)\h v (iD a )(iD f ,)(iD p )'js'l5h v \H 00 (v)) = \d H p 2 ie uol p 5 v u v IJi . (14) 



^They are related to the matrix elements p^ and p\ s introduced in Ref. |Tl[] by p\ = p%,p2 



3PLS- 
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The expectation value of any bilinear operator with three derivatives is expressible in terms 
of pi and p 2 : 



(HMlKTiiD^iiD^iiD^KlH^v)} = 
\pi {g a p - v a vp) v^Tr [P+Y] - j^d H p 2 ie va p & v v v^Tr P + ^ s ^ b P + T 



(15) 



where P + = ~(1 + and T is any four- by-four matrix. 

After a rather lengthy calculation we obtain the contributions from local dimension-six 
operators to the form factors: 



(3) 



+ 



Pi + 3p 2 _J_ 
12A m 2 b 2A% 
2(pi + 3p 2 ) 



Pi - P2 + 



(pi + 3p 2 ) (q 2 - q-v 2 - m| + m b q-v) 



3m 2 b 



SA d m b 



(m b - q-v) ( 



2 2 

q — q-v 



4p 



3A^ 



\{m b - q-v) 2 ( 



2 2\ 

q — q-v I 



(16) 



T, 



(3) _ Pi + 3p 2 



+ 



6A mj; 3Ag 
2 



4pi + 6p 2 - (pi + 3p 2 ; 



3A^ 



(4pi + 6p 2 )(m 6 - q-v)q-v - 3p 2 q" 



3A% 



(m b - q-v)(q 2 - q-v 2 ), 



(17) 



(3) 



Pi + 3p 2 



6m 2 Aq ^ 1 



2(m;, — g-v) 
3Al 



q-v 

(Pi + 3p 2 ) 3p 2 

m b 



- ^(m b - q-v)(q 2 - q-v 2 ), (18) 



(3) 



pi + 3p 2 2p 2 4(pi + 3p 2 



3m 2 A 2 



A 3 



+ 



3m b Al 



(m b - g-u) 



(19) 



(3) 



Pi + 3p 2 
6m 2 Aq 



q-v 



2(pi + 3p 2 ) 2p 2 m fe 4pi 2 2 

g-u (m 6 - q-v) + + ^-rj{m b - q-v)(q - q-v ). 



3m b Al 



A 3 



3A* 



(20) 



Substituting the imaginary part of these form factors into eq. (||) we obtain the corrections 
to the triple differential decay rate. Interesting quantities are the charged lepton spectrum 
and the hadronic spectrum. The former is obtained by taking the imaginary part of form 
factors T^\i = 1,2,3 and integrating eq. (H) over q 2 and E v . Using the rescaled lepton 



energy y = 2Ep/m b we find the l/m b correction to the lepton spectrum 



l-r-y) 



(3px + r 2 pi + 9r 2 p 2 - 3r 3 p 2 ) + -pi2/(3 - 2r) - 8/V (21) 



2. , 3 2 (4 Pl + 3r 2 Pl + 9r 2 p 2 ) 2r (8 Pl + 9r Pl + 27rp 2 

- 3 (pi + 3p 2 ) y 



+ 



1-y 3(1 -yf 

2r (8pi + 17rpi + 4r 2 pi - 9rp 2 + 12r 2 p 2 ) ^_ 2r 2 (3pi + 4rpi + 9p 2 + 12rp 2 ) 

3(1 -yf 

8r 2 (px + 3rpx + 3rp 2 ) 40r 3 pi 



(1 " y) E 



3(1 - 



5(1 - r - y) 



2(l-r) 4 (l + r) 2 pi 
3^2 



where r = (m c /mb) 2 . Note the contribution from the (^-function at the endpoint of the lepton 
spectrum. For b — > u transition such singular terms in the lepton spectrum appear already 
at order 1/m 2 , but for b — > c they do not appear until order 1/m 3 . This is easily explained 
if one recalls that the most singular contributions to the lepton spectrum at a given order 
1/m^ can be obtained from the spectator model result by the "averaging" procedure of 
Ref. 0, which involves differentiating n times with respect to y. For a massless final state 
quark the spectator model spectrum has the form f(y)8(l — y) with /(l) ^ 0, and thus 
differentiation produces the n — 1-st derivative of the ^-function 5( n_1 )(l —y). For a massive 
quark in the final state the spectator model spectrum and its first derivative vanish at the 
end point y = 1 — r. Hence at order l/m% the most singular contribution is proportional to 
*M(i- y - r ). 

To obtain the contribution from local dimension-six operators to the hadronic spectrum 
we integrate eq. (||) over E v and express the result in terms of rescaled hadronic variables 
Eq = (nib — q ■ v)/rrib and §o = (ml — 2q ■ v + g 2 )/m 2 : 

dr® sr 



ds dE 3m 3 , 



Q(E - Js )Q(l + s - 2E )^E 2 - s x 



(22) 



|(pi + 3p 2 )(-3 + 6E + 2E 2 - 5s )5(s - r) - 2[9{ Pl - p 2 ) 
+ 6(pi - p 2 )s + 3(pi + 3p 2 )So - (3(7pi - 3p 2 ) + ll(pi + 3p 2 )s )E 
+ 3((3 Pl + 5p 2 ) - ( Pl + 3p 2 )s )E 2 + 8(px + 3p 2 )E$]5'(s - r) 
- 4(E 2 - s ) [3(1 + s )p 2 - (1 - 3s )(pi + 3p 2 )£ - 2( Pl + Qp 2 )E 2 ]5"(s - r) 



+ -E (E 2 - s )pi 2S - 3(1 + s )E + AE 2 Q 6"'(s - r 



The correction to the total rate is given by integrating eq. (|21|) or eq. ([22|) over the remaining 
variables: 

p(3) = _0_ ^ _ gg r + 24r 2 _ gr 3 _ 5r 4 + 4g j og r + 36r 2 j og r ) ^3) 

+p 2 (27 - 72r + 216r 2 - 216r 3 + 45r 4 + 108r 2 log r)]. 

The part of eq. ( p3|) that diverges logarithmically as r — > agrees with the corresponding 
expression in Ref. There is nothing wrong with the logarithmic divergence, since 



our calculation is valid only for the charm mass significantly larger than Aqcd- It is the 
latter condition that allowed us to discard the diagram in Fig. lb. For a discussion of the 
corrections to the total semileptonic decay rate from dimension-six operators with a light 
quark in the final state see Ref. |T3| . 



III. EXPANSION OF THE STATES 

Above we have computed the \jm\ corrections to the inclusive differential B decay rate 
from the local dimension-six operators in the OPE. However, there are other sources of 
1 / ml corrections. At order 1 / m 2 the OPE yields the decay rate in terms of the two matrix 
elements 

(B(v)\h v (iD ± ) 2 h v \B(v)), (24) 
l(B{v)\h v l <v G>» h v \B{v)), 

where \B(v)) is the physical 5-meson state, rather than the state of the effective theory in 
the infinite mass limit [B^v)). Thus these matrix elements are mass-dependent. At order 
1/m 2 this distinction is irrelevant, but at higher orders this mass dependence has to be taken 
into account explicitly. We express the physical states through the states in the infinite mass 
limit of HQET using the Gell-Mann and Low theorem (see e.g., Ref. Jl4[). This theorem 



implies that, to first order in l/m?,, \B(v)) is given by 

r° . , s 1 



B(v)) 



l + i [d 3 x [ dtCi{x) - -(B^v) \i f d 3 x f dtCi(x)\ Booiv)) 

J J —CO V J J —OO 



Boo(v)), 

(25) 



where V is the normalization volume and 

Ci = ^-hv (iD ± ) 2 h v + ^-h v | <v GT h v . (26) 

Utilizing eq. (|25|), one can easily expand the matrix elements in eq. ( |2^|) to order l/m\. It 
is convenient to introduce the following notation: 



{H^KiiD^Ki f d 3 x f° dt£ I {x)\H 00 (v))+h.c. = Tl + dHT \ (27) 
(H^hJ-a^G^Ki f d 3 x f° dtC I (x)\H 0O (v)) + h.c. = - + — . 



We then find 



T 4- ^7" 

K (^ ± ) 2 ^| B{v)) = Ax + , (28) 

i |^ | cv fc„| S(«)> = A 2 + T3 Q +37 " 4 . 
3 2 3m fe 

Thus these order 1/mf corrections to the inclusive B —>■ XJLv decay rate are parametrized 
by the matrix elements T\ — of four nonlocal operators^ 

This class of 1/mjj corrections can be included in any quantity known at order l/mj^ by 
using eq. ( P8"D to evaluate the matrix elements of the dimension-five operators. In particular 
the corrections to the form factors and the differential rates in Ref. || can be obtained in 
this way. 



IV. APPLICATIONS 

One important application of our results is to study the influence of 1/mf corrections on 
the extraction of the HQET matrix elements A, Ai using the methods of Refs. |8],|9|]. 

In order to compare quantities obtained from an expansion in the inverse quark mass 
with experiments it is necessary to express the quark masses m c and mj, through the physical 



^These matrix elements are related to those introduced in Ref. [11] as T\ = p 3 T7V ,T2 = ^p^.Q,7s 



G- 
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meson masses mg and mo and the HQET matrix elements. Some details of this calculation 
are given in the appendix. To order 1 / m| we find the following relation 

- Ai + d H \ 2 (m Q ) pi + d H p 2 T 1 + % + d H (T 2 + T 4 ) 
'""= m « +A 2^ + i^J • < 29 > 

where m B is the hadron mass and rug is the heavy quark mass. The differential and total 
decay rates are functions of the ratio of quark masses which can be expressed in terms of 
the spin averaged meson masses 

m E= ^D_^_/ 1 _mn\^rmn_mn j » (30) 
m b m B m B \ m B J 2m B \m D vtib "" "' 
A 3 / rn B) \ AAi / m B 



A 2 


(- 










m B J 




m%\ 


Pi 




i ( m 2 B 


m 2 D ) 






\m 2 D 



_< ,l-= il +^- l + =^-3 = ^ + 
V mfi/ irrfg y m B m B 

where m D and m B are defined as rfiMeson = ( m P + 3my)/4. 

The familiar relation of the HQET matrix element A2 to the mass splitting between B 
and B* mesons also needs to be extended to include the \jm\ contributions. Using eq. ( |29|) 
to express the quark mass through the meson mass and A, we find 

n(m Q )\ 2 (m b ) , A x p 2 T 2 + % 

m H * - m H = Am H = 2 (H 5- H 5 — , (31) 

rriH v rr%H n^H m H 

where ^{mq) = (a s (mQ) /a s (mb)) 3 ^ takes account of the scale dependence of A2. We can 
use the B — B* and D — D* mass splitting to extract the numerical value of some of the 
HQET matrix elements: 

Am B m B — Am B 'm 2 D 



A 2 (m&) 



2(m B - K,(m c )m D ) 



K(m c )m 2 B Am B (m D + A) - m 2 D Am D (m B + A) 

p 2 — 1 2 — i4 — = p — ; tT . (OZ) 

mg + A — K[m c )[m B + A) 

In order to extract Ai and A from the experimentally measured lepton energy spectrum 
in the B — > X c £v decay it is convenient to introduce the quantities |§ 

r> _ J 1.5 GeV dEe D _ 7l.7GeV dE^ / q o \ 

Rl - — dT~ ' R2 - 7 dT~~ ' (33) 

■ dEi I ttt~ dE e 



1.5 GeV dEe J 1.5 GeV d-E 
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where Ei is the lepton energy and dT/dEe is the complete electron energy spectrum, which we 
obtain by combining the results of Ref. H with our results. In the energy spectrum at order 
1/ml, taken from Ref. [[J, the matrix elements of dimension-five operators are evaluated 
according to eqs. fl2"8"|). The resulting expression is combined with the contribution from local 
dimension-six operators eq. (0). Expressing all quark masses through the meson masses 
and using their measured values, we obtain expressions for i?i,i?2 in terms of the HQET 
matrix elements. Combining these with perturbative corrections and other contributions 
(see Ref. §) we find 

A A 2 A A A 3 AA 

iMGeV] = 1.8059 - 0.309 — - 0.35 ^ - 2.32 - 3.96 ==§- - 0.4— 5- - 5.7=^ 

m B m B 171 B 171 B m B m B 

-6.8===^ - 7.7==-, - 1.3=5- " 3.2=| - 4.5^. - 3.1=5- - 4.0=5, 

m B m B m B m B m B m B m B 



— (0.035 + 0.07 ==M + 

7r V m B J 



V, 



ub 



cb 



2 



1.33 - 10.3 — | - (0.0041 - 0.004 — 
m B V m B/ 



+ (0.0062 + 0.002 ==M , (34) 

V ™>b) 

A A 2 A A A 3 Aa 

R 2 = 0.6581 - 0.315 — - 0.68 =^ - 1.65 — - 4.94 =§- - 1.5 - 7.1 — ^ 

m B m B 171 B 171 B m B 171 B 

AAo Pi Pi % 7o 7q Ta 

-17.5 =^ - 1.8 ±5r + 2.3 ±5r - 2.9=^- - 1.5=V - 4.0=-!- - 4.9=^ 

m B m B m B m B 



— (0.039 + 0.18==^- 

7T V 



u6 



V, 



cb 



2 



0.87 - 3.8 — ) - (0.0073 + 0.005 ===- 

m B l \ m B , 



+ (0.0021 + 0.003 ==M , (35) 

where the first two lines contain the nonperturbative corrections to order l/fn B . The other 
terms are in order: the perturbative a s corrections, the contribution from B — > X u lv decays, 
electroweak corrections, and finally a boost, since the 5-mesons do not decay from rest. 
This is to be compared with the experimental values R exv = 1.7831GeV, R^ xp = 0.6159. 
Neglecting the l/m B corrections but including statistical errors the values A = 0.39 ± 
O.llGeV and Ai = — 0.19+0. lOGeV 2 were found in Ref. 0. In order to take the uncertainties 
from the higher order matrix elements into account, we equate the expressions for i? 12 to 
the experimental values using \V u b/V c t\ = 0.08, a s = 0.22 and eqs. ( |3"2"D to eliminate A 2 and 

12 
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FIG. 2. Impact of l/mf corrections on the extraction of A, Ai- Shaded region: Higher order 
matrix elements estimated by dimensional analysis. Cross-hatched region: p\ = 0.13GeV 3 , p 2 = 0. 
Cross and ellipse show the values of A, Ai extracted without corrections but including the 

experimental statistical error. 

p 2 . This yields the extracted values of A, Ai in the form 



A — f/ i (Rl xp ,Rl xp ,pi,T 1 ,T 2 ,T 3 ,T 4 ), Ai — f Xl (Rl xp , Rl* 1 ' , pi, Tx, T 2 , %, T 4 ) 



(36) 



Dimensional analysis suggests that the higher order matrix elements are all of order 
Aq CD , which can be used to make a quantitative estimate of the uncertainties in the ex- 
traction of A, Ai. We vary the magnitude of pi,7{ — 1i in eqs. ( |56] ) independently in the 
range — (0.5GeV) 3 , taking p\ to be positive, as indicated by the vacuum saturation ap- 
proximation, but making no assumption about the sign of the other matrix elements. Using 
the central values for R ex 2 v we find that A, Ai can lie inside the shaded region in Fig. 2. For 
comparison we also display the values of A, Ai extracted in Ref. |§ together with the ellipse 
showing the size of the statistical error of the experimental data. Clearly the theoretical 
uncertainties dominate the accuracy to which A, Ai can be extracted. 

The situation can be improved only if we have some independent information on some or 
all of the higher dimension matrix elements. This requires either more experimental input or 
theoretical estimates of these matrix elements. p\ can be estimated in the vacuum saturation 
approximation [ ]T^|n||n] , |r^|r7| || , p 1 = (2na s / '^mg/^. The numerical value obtained this 
way is rather uncertain. Taking a s = 0.5 and fs = 270MeV for purposes of illustration, we 
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find pi ~ O.13GeV 3 .0 No similar estimates exist for the other dimension-six matrix elements. 
P2 vanishes in any non relativistic potential model, which may be taken as an indication that 
it is small relative to the other matrix elements. No estimates that go beyond dimensional 
analysis are available for the time ordered products. 

The cross hatched region in Fig. 2 shows the range of A, Ai one obtains from setting 
pi = 0.13GeV 3 and P2 = and varying the magnitude of the other matrix elements in the 
range — (0.5GeV) 3 . The previously extracted values of A, Ai are not excluded by this choice 

of Pl,2- 

This method of extracting A, Ai is especially sensitive to higher order corrections since 
the constraints obtained form R\ and R2 give almost parallel bands in the A — Ai plane. 
Thus small uncertainties in the theoretical expressions for Ri 2 result in large uncertainties 



in the extracted values of A, Ai. The same applies to the very similar analysis in Ref. |L8 
The rare decay B — > X s 7 provides a way to extract a vertical band in the A — Ai plane, but 
at present the experimental data does not allow a quantitative analysis |[19|| . Furthermore, 



as discussed in the introduction, it is not clear when HQET matrix elements extracted from 
different observables can be compared meaningfully 

The second method for extracting information on A, Ai [[| was used to exclude some 
regions in the A — Ai plane. The first and second moments of the invariant mass spectrum of 
the hadrons in the final state of the inclusive decay B — > X c iv turn out to give independent 
constraints on A, Ai. Their definition involves the total decay rate at order l/m 3 B . It can 
be obtained by combining the total rate at order 1 / m| from Ref. || with the contributions 
from local dimension-six operators eq. (|23|) and using eqs. (p8|). Finally eqs. ( j2"9|) and (|30D 
are used to eliminate the quark masses. Using the measured values for the meson masses 
and neglecting perturbative corrections we find to third order in l/m^: 



= A value for pi can also be obtained from small velocity sum rules [17] but this estimate suffers 
from large uncertainties as well. 
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\V cb \ 2 G 2 F m% 



1.21 



1927T 3 

AAi 



A A 2 A A A 3 

0.3689 - 0.6080— - 0.349^ - 1.175— J- - 2.757^- - 0.11=^ (37) 
m B m B m B m B m B 



A \ T" T" T" T 

2.95^ - 2.27-^- + 2.76-^- - 2.73^- + 0.55^- - 3.84^- - 2.76^- 



m 



l b ^1>B llv B llfj B llfj B llv B 11 V B ltL B 

Since none of the coefficients of the higher order matrix elements turn out to be abnormally 
large, dimensional analysis indicates that the l/m B corrections to the total rate should not 
exceed 2%. 

The hadronic moments are defined as 

dr 



mi 



m 



mi 



mi 



mi 



((s H - m 2 D ) n ) = y J ds H dE H (s H - m 



—2 \n 
D) 



(38) 



ds H dE H 

where sh = rn 2 B — 2m sv-q + q 2 and E B = ms — v-q are the hadronic analogs of so , E defined 
in Sect. II. Using the relation between quark and hadron masses one can relate sh, E b to 
So, E Q and thus compute the moments using the expressions given in Ref. |§ together with 
eq. (|2~2"1) and the usual substitution eqs. (^). We find to order l/m B : 



(s H - m 2 D ) = m B 



a q A / a« 

0.051— + 0.23— 1 + 0.43— 

7r m B \ IX 



0.26 =;r (A 2 + 3.9A 1 - 1.2A 2 ) 



m 



B 



+0.33 =;r (A 3 + 6.6AA1 - 1.7AA 2 + IJQpx + 3.5p 2 
m B 



+5.0T! + 2.5T 2 + 4.6T 3 + 1.3T 4 ) 



{{sB-rn'nY) 



m 



B 



0.0053— - 0.038—— + 0.065^- (A 2 - 2.1Ai) 



.OCs 
7T 



m B it 



m'i 



1 



+0.14 = ^(A 3 + 2.2AA X + 2.2AA 2 - 6.O01 + 1.7p 2 - l.OTi - 2.9T 2 ) 
m B 



(39) 



(40) 



where perturbative a s corrections have been included. Rather that repeating the analysis 
presented in Ref. , we use these expressions to predict the values of the hadronic moments 
using the HQET matrix elements extracted from the lepton energy spectrum. The main 
reason for doing this is that the experimental measurement of the necessary branching 
fractions is not very precise. In particular ALEPH and CLEO quote only an upper bound 
for Br(B — > D\iv) pO| , |2T| . We extract an upper bound on this branching fraction from 
the theoretical prediction of the hadronic moments. A lower bound for the first hadronic 
moment is given by |§ 
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(s H - m D ) > a 



(2.450GeV) 2 - (1.975GeV) s 



+ b 



(2.010GeV) 2 - (1.975GeV) s 



+c 



;i.869GeV) 2 - (1.975GeV) s 



(41) 



where a, b, and c are the semileptonic branching fractions to D**, D*, and D relative to the 
total semileptonic branching fraction. Using the measured ratio 0.41:0.59 for the decays to 
D and D*, we can write b, c as functions of the branching fraction a for D**: 



6 = 0.59(1 -a) 



0.41(1 - a) 



(42) 



We take a + b + c = 1, which is appropriate because we need only a lower bound on the 
hadronic moment. It is also implicitly assumed that the nonresonant semileptonic branching 
fraction below the D** mass is negligible. Similarly, for the second hadronic moment we take 



{(s H - m D ) ) > a 



(2.450GeV) 2 - (1.975GeV) 5 



(43) 



where small contributions from the ground state mesons D, D* have been neglected. We 
obtain theoretical predictions for the hadronic moments by substituting values of p±, T\ — 71 
and the corresponding values of A, Ai extracted from the lepton spectrum into eqs. (|39|) . As 
before, we allow the magnitudes of p\ and T\ — to vary in the range — (0.5GeV) 3 with 
pi being positive. Imposing the constraint that the largest values of the hadronic moments 



obtained from this procedure be larger than the lower bounds eqs. ([|l]),([I^) we find the 
upper bound on the D** branching fraction 

a < 0.23. (44) 

This value is compatible with the experimentally measured values from ALEPH |2(| 
(Br(B -> D x iv) = 0.069 ± 0.015, Br(B -> D*£v) < 0.11) and from CLEO ^(Br(B -> 
D-JLv) = 0.046 ± 0.013, Br(B -> D* 2 iv) < 0.11). It is also marginally consistent with the 
OPAL result a = 0.34 ±0.07 ||22|| . Unless the matrix elements of dimension-six operators are 
even bigger than we have assumed, this implies that the branching fraction a = 0.27 used 
in II is inconsistent with the values of A, Ai extracted from the lepton spectrum. 
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V. CONCLUSIONS 



We have calculated the l/m\ contributions to various observables in the semileptonic 
decay B — > X c iv. They are parametrized by the expectation values of two local and four 
nonlocal dimension-six operators. While the total rate is rather insensitive to the higher 
order corrections (1-2%), the values of A, Ai extracted from the lepton spectrum can be 
affected substantially. The theoretical uncertainties in the values of A, Ai are far larger than 
the statistical errors of the experimental measurements if the values of the higher order 
matrix elements are estimated using dimensional analysis. While one linear combination of 
A and Ai is still reasonably well constrained, it is not possible to extract individual values for 
A and Ai from the lepton spectrum only. The situation can be improved only if additional 
information on the size of the dimension-six matrix elements is used. Unfortunately no 
theoretical estimates are available for any of these matrix elements except p\. The latter 
can be estimated in the vacuum saturation approximation, albeit with large uncertainties. 
Alternatively one can use additional experimental input, e.g., from B — > X s 7 decays, to 



further constrain A and Ai [pLSf] . 

The values of A, Ai extracted from the lepton spectrum can be used to make theoret- 
ical predictions for the moments of the hadronic invariant mass spectrum. This amounts 
to expressing one observable in terms of other observables, a procedure that makes sense 
only if the perturbative series for this expression is reasonably well behaved. In order to 
determine whether this is the case it is necessary to know at least the next-to-leading or- 
der a s corrections to all observables involved. Since they have not been computed for the 
lepton spectrum, there is at present no way we can check if predictions for the hadronic 
moments in terms of the HQET matrix elements extracted from the lepton spectrum satisfy 
this criterion. 

Setting these considerations aside, we can predict the values of the hadronic moments 
in terms of A, Ai extracted from the lepton spectrum. The lower bounds for these moments 
depend on the branching fraction to D**, which is not well known experimentally. By 
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demanding that not the whole range of predicted values of the hadronic moments be excluded 
by the lower bounds, we find an upper bound of 23% on the branching fraction to D**, if the 
higher order matrix elements are estimated by dimensional analysis. This value is consistent 
with the ALEPH and CLEO measurements. 
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APPENDIX A: THE MASS FORMULA 

For comparison with experiments it is necessary to express the pole quark masses m c 
and mi, in terms of HQET matrix elements and physical observables, e.g., the spin averaged 
meson masses m B and m D , where m Meson = {m P + 3my)/4. For this purpose one needs to 
know how quark masses are related to hadron masses at order l/m\. Our starting point is 
the identity 



m g = - x ; r "7 wl T ,\r +h.c., (Ai) 



V(H OQ {v)\ 


h\h(v)) 


2 {H 00 {v)\ 


H(v)) 



where V is the normalization volume and 7i is the full Hamiltonian density including light 
degrees of freedom^. This equation holds in the rest frame of the hadron. Then we split H, 
into the leading term and the terms suppressed by powers of l/m b , H = 7Y + 'Hi, and use 
the fact that | H^v)) is an eigenstate of / d 3 xHo with eigenvalue m b + A. The use of the 
Foldy-Wouthuysen-transformed fields eq. ( [T0| ) ensures that there is no implicit dependence 
on nib in h v . Also, there are no time-derivatives in the HQET Lagrangian beyond leading 



**If one starts from a similar identity with eigenstates of TL on both sides of the matrix element, 
one obtains the same result after a somewhat more cumbersome calculation. 
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order, as can be seen e.g. from eq. (82) of Ref. Therefore we have Tii = —C\. Using 



the Gell-Mann and Low theorem, the general expression for the hadron mass reads 



x V 
m b + A - — 





C a Texp (ifcPx^dtC^x)] 


) 1 HM) 




Texp (i J d 3 x f®^ dtC\{x)\ 





+ h.c. 



Expanding eq. (|A2|) to order 1/mf we obtain the mass formula: 

m H = m b + A - (H^v) \d + C n \ H^v)) 

Uh^v) \Cji [ d*x [° dtCj(x)\ H^v)) + h.c. 



where []T0 

c n = - 



4m? 



1 1 1 - 

h v ip±(iv ■ D)ipj_ h v + ■^^hy (ip±) 2 (iv ■ D) h v + - — ^h v (iv ■ D){ii 



■smt 



imi 



(A2) 



(A3) 



\) 2 h v , 



(A4) 



and Cj is given in eq. fl26|) . Eq. (|A3|) contains expectation values of both local and nonlocal 
operators. The local part can be evaluated in terms of the matrix elements Ai, A2, P\ and p 2 , 
while the nonlocal matrix elements can be expressed through 71 — 74 defined in eqs. 
In terms of these matrix elements the meson mass is given by 



-j X 1 + d H X 2 pi + d H p 2 71 + 7j + d H (T 2 + %) 
m H = m b + A h 



2m b 



4ml 



Ami 



(A5) 



in agreement with Refs. [10,11 . 
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